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We  study  a  two-state  statistical  process  with  a  non-Poisson  distribution  of  sojourn  times.  In  accordance  with 
earlier  work,  we  find  that  this  process  is  characterized  by  aging  and  we  study  three  different  ways  to  define  the 
correlation  function  of  arbitrary  age  of  the  corresponding  dichotomous  fluctuation.  These  three  methods  yield 
exact  expressions,  thus  coinciding  with  the  recent  result  by  Godreche  and  Luck  [J.  Stat.  Phys.  104,  489 
(2001)].  Actually,  non-Poisson  statistics  yields  infinite  memory  at  the  probability  level,  thereby  breaking  any 
form  of  Markovian  approximation,  including  the  one  adopted  herein,  to  find  an  approximated  analytical 
formula.  For  this  reason,  we  check  the  accuracy  of  this  approximated  formula  by  comparing  it  with  the 
numerical  treatment  of  the  second  of  the  three  exact  expressions.  We  find  that,  although  not  exact,  a  simple 
analytical  expression  for  the  correlation  function  of  arbitrary  age  is  very  accurate.  We  establish  a  connection 
between  the  correlation  function  and  a  generalized  master  equation  of  the  same  age.  Thus  this  formalism, 
related  to  models  used  in  glassy  materials,  allows  us  to  illustrate  an  approach  to  the  statistical  treatment  of 
blinking  quantum  dots,  bypassing  the  limitations  of  the  conventional  Liouville  treatment. 
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I.  INTRODUCTION 

The  phenomenon  of  aging  has  been  known  for  a  long  time 
to  be  a  property  of  spin  glasses  and  polymers  [1].  Part  of  the 
reason  for  the  recent  interest  in  the  aging  phenomenon  has  to 
do  with  the  predicted  breakdown  of  certain  fundamental  as¬ 
sumptions  made  in  equilibrium  statistical  mechanics  when 
applied  to  strongly  disordered  systems.  For  example,  the  On- 
sager  principle  [2],  which  is  the  relaxation  of  a  perturbed 
system  back  to  its  equilibrium  state  described  by  an  unper¬ 
turbed  autocorrelation  function,  is  violated  in  anomalous  dif¬ 
fusion  leading  to  an  anomalous  relaxation.  More  recent  pa¬ 
pers  on  this  phenomenon  are  devoted  to  studying  aging  in 
diffusion  processes  occuring  in  of-dimensional  lattices  [3],  in 
low  dimensional  environments  [4]  and  in  the  quantum  dy¬ 
namics  of  dissipative  free  particles  [5],  Most  recently  [6,7] 
there  has  been  some  interest  in  the  manifestation  of  aging  in 
processes  described  by  means  of  the  continuous  time  random 
walk  (CTRW)  formalism  [8]. 

The  recent  work  of  Ref.  [9]  establishes  that  the  intermit¬ 
tent  radiation-induced  fluorescence  of  colloidal  semiconduc¬ 
tor  quantum  dots  also  yields  aging.  These  new  materials, 
under  the  influence  of  a  radiation  field  generate  intermittent 
fluorescence  suggesting  that  the  system  jumps  back  and  forth 
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from  a  “light  on”  to  a  “light  off'  state  [10],  with  a  non- 
Poisson  distribution  of  times  of  sojourn  in  these  two  states. 
This  is  a  quantum  mechanical  process  that,  in  principle, 
should  be  described  using  a  two-state  master  equation.  How¬ 
ever,  the  important  work  of  Jung  et  al.  [10]  shows  that  the 
adoption  of  the  ordinary  Liouville  approach  to  the  density 
time  evolution  encounters  serious  difficulties,  caused  by  the 
presence  of  aging.  This  difficulty  is  confirmed  by  the  more 
recent  work  of  Aquino  et  al.  [11],  who  have  evaluated  the 
emission  and  absorption  spectrum  of  this  intermittent  system 
using  the  CTRW  formalism.  For  this  reason,  it  is  important 
to  address  the  problem  of  the  foundation  of  the  master  equa¬ 
tion  from  the  new  perspective  of  trajectories  (CTRW)  rather 
than  from  that  of  densities,  as  in  the  fundamental  work  of 
Zwanzig  [12].  The  conflict  between  these  two  pictures,  no¬ 
ticed  for  the  first  time  by  Bologna  et  al.  [13],  is  far  from 
being  fully  understood,  and  it  is  not  yet  clear  if  the  contro¬ 
versy  is  due  to  the  inadequacy  of  the  known  forms  of  non- 
Markov  master  equations,  derived  from  the  Liouville-density 
method,  or  to  an  even  more  fundamental  reason,  namely,  the 
breakdown  of  the  generally  accepted  equivalence  between 
the  density  and  trajectory  pictures.  To  shed  light  into  this 
important  issue,  herein  we  plan  to  reexamine  the  issue  of 
aging,  widely  discussed  in  the  field  of  glassy  materials,  and 
the  important  mathematical  treatment  of  this  property  made 
by  Godreche  and  Luck  [14]  as  well,  from  the  perspective  of 
the  generalized  master  equation  (GME). 

A  natural  way  of  formally  expressing  aging  is  through  the 
correlation  function  of  the  stochastic  variable  of  interest,  f, 
which  is  expressed  as  an  ensemble  average  indicated  by  the 
brackets 
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<^l)^2)>,a 

<£2> 


=  <S>f\h-t2).  (1) 


As  we  shall  see  throughout  this  paper,  the  parameter  ta  de¬ 
notes  the  time  distance  between  the  preparation  and  the  ob¬ 
servation  process.  The  averaging  brackets  carry  a  subscript 
stressing  that  this  is  a  f„-old  property,  rather  than  the  tradi¬ 
tional  aged  property,  only  depending  on  |fi-f2|.  The  time 
difference  indicates  the  stationarity  of  the  underlying  pro¬ 
cess,  but  the  ta  superscript  denotes  a  dependence  on  the  age 
of  the  process  as  well.  In  the  case  of  a  dichotomous  variable 
obeying  renewal  theory,  the  exact  expression  for  the  Laplace 
transform  of  this  age-dependent  correlation  function  was 
found  by  Godreche  and  Luck  [14].  Their  exact  result  was 
recently  recovered  by  Margolin  and  Barkai  [15]  as  a  special 
case  of  a  more  general  expression,  since  the  latter  authors  do 
not  require  the  condition  that  the  two  states  of  the  dichoto¬ 
mous  variable  i  have  the  same  waiting-time  distribution. 
However,  herein  we  make  the  same  assumption  as  did  Go¬ 
dreche  and  Luck  [14],  and  derive  their  result  along  the  lines 
of  the  recent  work  of  Ref.  [7].  Allegrini  et  al.  [7]  noticed  that 
in  the  non-Poisson  case  the  well  known  GME  of  Kenkre, 
Montroll,  and  Shlesinger  [16]  becomes  incompatible  with 
the  Onsager  principle  [2]  and  found  a  way  to  make  the  GME 
compatible  with  the  aged  condition.  However,  they  [7]  left 
unsolved  the  problem  of  deriving  a  GME  of  arbitrary  age, 
which  would  be  equivalent  to  making  the  Onsager  principle 
compatible  with  an  incompletely  aged  system.  Herein  we 
solve  this  problem  and  establish  that  this  solution  leads  to  the 
exact  expression  for  the  correlation  function  of  arbitrary  age 
found  by  Godreche  and  Luck  [14]. 

The  GME  derived  herein  does  not  have  the  same  origin  as 
those  widely  discussed  in  the  literature  [12].  The  Zwanzig 
GME,  in  fact,  one  of  the  most  popular  master  equations,  is 
derived  from  a  first  principles  procedure,  starting  from  the 
statistical  Liouville  equation  of  the  whole  universe.  The 
Zwanzig  GME  is  a  projection  of  this  universal  Liouville 
equation  onto  the  Hilbert  space  of  the  system  of  interest.  The 
approach  used  to  derive  the  GME  herein  is  based  on  the 
experimental  observation  of  non-Poisson  dichotomous  sig¬ 
nals.  Examples  of  such  signals  include  those  produced  by 
ionic  channels  [17]  and  by  blinking  quantum  dots  [10].  We 
build  up  a  GME  that  is  compatible  with  the  experimentally 
determined  non-Poisson  nature  of  these  processes,  assuming 
the  applicability  of  renewal  theory.  We  leave  open  the  ques¬ 
tion  as  to  the  source  of  randomness,  but  note  in  passing  that 
the  fluctuations  in  the  system  variable  are  generated  by  the 
environment.  A  simplified  dynamical  model  has  reproduced 
the  essential  statistical  properties  of  the  system  of  interest,  in 
the  case  of  infinitely  aged  GME  [18].  In  the  case  of  arbitrary 
age,  not  only  a  first  principle  derivation  is  missing,  a  condi¬ 
tion  shared  by  the  infinitely  aged  case:  there  is,  to  date,  not 
even  a  heuristic  derivation  of  the  corresponding  GME. 
Herein  we  provide  a  heuristic  derivation  of  a  GME  with 
arbitrary  age,  one  based  on  an  exact  expression  for  the  cor¬ 
responding  distribution  of  sojourn  times,  and  we  use  this 
result  to  derive  an  exact  expression  for  the  correlation  func¬ 
tion. 


In  concluding  the  Introduction,  we  point  out  that  we  re¬ 
work  the  problem  from  a  number  of  perspectives.  We  do  this 
because  each  approach  provides  separate  and  distinct  in¬ 
sights  into  the  phenomenon  described  by  a  GME  of  arbitrary 
age.  The  first  perspective  adopted  in  this  paper  focuses  on 
the  derivation  of  a  GME,  based  on  the  experimental  obser¬ 
vation  of  the  time  evolution  of  a  trajectory,  characterized  by 
rare  jumps  from  one  state  to  another.  The  second  perspective 
uses  a  Liouville-like  approach  to  the  time  evolution  of  the 
variable  of  interest  and  its  environment.  No  attempt  is  made 
to  derive  the  former  dynamic  picture  (trajectory)  from  the 
latter  (probability  density).  The  probability  density  perspec¬ 
tive  is  characterized  by  infinite  memory,  yet,  the  statistical 
process  under  study  is  generated  by  random  critical  events, 
whose  occurrence  erases  the  memory  of  earlier  events. 
Therefore,  we  think  it  is  prudent  to  examine,  yet  again,  the 
same  process  from  a  third  perspective,  one  based  on  the  di¬ 
rect  observation  of  the  sequence  of  rare  random  events. 

We  find  that  the  use  of  three  different  perspectives  is  fruit¬ 
ful  in  shedding  light  on  the  recent  observation  made  by 
Sokolov,  Blumen,  and  Klafter  [19].  Their  results,  in  our 
opinion,  imply  the  breakdown  of  certain  well-established  no¬ 
tions  of  equilibrium  statistical  mechanics,  such  as  linear  re¬ 
sponse  theory  as  a  prescription  for  predicting  the  effect  of 
externally  perturbing  a  system  out  of  equilibrium. 

The  outline  of  the  paper  is  as  follows.  In  Sec.  II  we  em¬ 
ploy  the  reduced  density  perspective,  namely,  we  build  up 
the  GME  of  arbitrary  age.  In  Sec.  Ill  we  study  the  time 
evolution  of  the  total  distribution  density.  In  Sec.  IV  we  use 
a  third  method,  based  on  the  probability  of  occurrence  of  the 
rare  random  events.  We  devote  Sec.  V  to  illustrate  the  physi¬ 
cal  connection  between  spin-glass  and  blinking  quantum  dot 
aging.  Finally,  in  Sec.  VI  we  make  some  concluding  re¬ 
marks. 


H.  THE  REDUCED  DENSITY  PERSPECTIVE 

This  section  is  devoted  to  the  application  of  the  reduced 
density  perspective  to  the  construction  of  a  GME.  In  the 
statistical  physics  literature  this  perspective  implies  a  con¬ 
traction  on  the  total  density  matrix  onto  a  prescribed  sub¬ 
space.  Consequently,  constructing  the  equation  of  motion  of 
the  total  density  matrix  should  be  the  first  step.  The  reason 
why  we  reverse  the  perspective,  thereby  confining  the  total 
density  treatment  to  the  next  section,  is  due  to  the  fact  that 
non-Poisson  statistics  has  the  effect  of  making  the  ordinary 
approach  extremely  difficult,  if  not  impossible,  as  discussed 
in  the  companion  papers  [18,20].  The  derivation  of  the  GME 
is  made  possible  by  expressing  the  higher-order  correlation 
functions  in  terms  of  the  second-order  correlation  function, 
via  simple  expressions  that  are  violated  by  the  non-Poisson 
statistics.  This  is  the  reason  why  we  derive  the  GME  from 
the  continuous  time  random  walk  (CTRW)  of  Montroll  and 
Weiss  [8,21],  rather  than  from  the  equation  of  motion  for  the 
total  distribution  density. 

A.  From  the  generalized  master  equation  to  the  age-dependent 
correlation  function 

The  reduced  density  matrix  perspective  is  based  on  the 
adoption  of  the  GME  formally  defined  by 
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yP (*)  =  -  f  ®(f-*')Kp(f')A\  (2) 

at  J0 


where  p(f)  is  the  //(-dimensional  population  vector  of  m  sites, 
so  that  its  ith  component,  /?,(/),  represents  the  probability  of 
finding  the  walker  at  time  t  in  the  ith  site,  states.  K  is  a 
transition  matrix  between  the  sites  and  <f>(f)  is  the  memory 
kernel.  We  study  the  case  where  the  fluctuating  variable  £ 
leading  to  this  GME  is  dichotomous.  This  means  that  the 
GME  has  only  two  states  and  the  matrix  K  reads 


The  waiting-time  distribution  in  either  of  the  two  states  is 
denoted  by  <p{t).  We  establish  the  connection  of  the  GME 
with  the  CTRW  [7,8]  by  relating  the  Laplace  transform  of 
the  memory  kernel  <!>(/),  given  by  fT>(/<),  to  the  Laplace  trans¬ 
form  of  the  waiting-time  distribution  ifAt),  given  by  <//(«),  as 
follows: 


<&(«)  = 


lit] Au) 

1  -  i/Au ) 


(4) 


the  first  laminar  region  is  extended  in  time  more  than  the 
time  interval  ta  between  preparation  and  observation.  The 
second  term  takes  into  account  the  cases  when  the  last  lami¬ 
nar  region  began  after  the  preparation  time,  after  a  sequence 
of  an  arbitrarily  large  number  of  earlier  laminar  regions,  the 
first  of  which,  of  course,  begins  at  t=-ta. 

In  the  Poisson  case,  because  of  its  unique  functional  form, 
there  is  no  dependence  of  i/f,  on  ta,  and  consequently  no 
aging.  In  the  non-Poisson  case,  on  the  contrary,  the  two 
waiting-time  distributions,  i/j,  and  i/A  t) ,  are  identical  only  if 
ta= 0.  In  this  case  both  CTRW  and  GME  correspond  to 
switching-on  the  observation  process  at  the  same  time  as  the 
preparation  process,  and  the  connection  between  the  two  pic¬ 
tures  is  given  by  Eq.  (4). 

Allegrini  et  al.  [7]  proved  that  the  GME  is  compatible 
with  an  infinitely  aged  CTRW,  provided  that  the  memory 
kernel  is  made  compatible  with  an  infinitely  aged  con¬ 
dition,  characterized  by  a  distribution  of  first  sojourn  times, 
which  is  infinitely  aged.  In  this  case  >°°  so  that 


where 


1  +  ifAu)  ~  2  $*,(«) 


(7) 


We  consider  two  distinct  processes,  that  of  preparation 
and  that  of  observation.  The  preparation  process  establishes 
the  initial  conditions  of  the  set  of  trajectories  under  study.  A 
trajectory  is  a  sequence  of  symbols  +  or  — ,  specifying 
whether  the  system  is  in  the  state  |  +)  or  |  — ).  We  call  the  time 
interval  with  the  system  either  entirely  in  the  state  |+)  or 
entirely  in  the  state  |— )  the  laminar  region.  With  Zumofen 
and  Klafter  [22],  we  assume  that  the  preparation  process, 
beginning  at  time  t=-ta,  insures  that  all  these  trajectories 
begin  with  the  system  at  the  onset  of  a  laminar  region,  either 
+  or  -. 

The  observation  process  begins  at  time  t=0^—ta.  The 
distribution  of  first  sojourn  times  is  denoted  by  <//,  (f)  and  is, 
in  general,  different  from  In  fact,  the  laminar  region 
corresponding  to  the  first  sojourn  time  might  have  begun 
earlier  than  t= 0.  The  exact  expression  for  this  time  distribu¬ 
tion  is 


<Af  ( t )  =  4At+t, 


00  rt 

I  a 

a)  +  2  I 

n=  1  J  0 


dy  i{Ay  + 1)  i[in{ta  ~  y) ,  (5) 


i  r 

dt'ifAt’)-  (8) 

\vJ , 

It  is  straigthforward  to  extend  the  calculations  of  Ref.  [7]  to 
the  case  of  an  arbitrarily  ffl-old  system,  so  that  Eq.  (7)  is 
replaced  with 


<M«)  = 


ui/jt  ( u ) 


1  +  <//(«)  -2  if/,  (n) 


(9) 


where  <//,  (u)  is  the  Laplace  transform  of  (//,  (f). 

Using  the  GME  with  the  ffl-old  memory  kernel,  we  define 
the  age-dependent  correlation  function  &g(t)  through  its 
Laplace  transform,  as  follows: 


= - - - . 

u  +  2<f>,  (m) 

a 

This  prescription  corresponds  to  setting 


(10) 


where  ipn{T)  denotes  the  probability  that  n  jumps  occur  dur¬ 
ing  the  time  interval  of  length  r,  the  last  of  which  occurs  at 
time  t=  t.  As  is  well  known,  for  a  renewal  process  the  wait¬ 
ing  times  for  successively  more  jumps  is  given  by  the  con¬ 
volution 


*/'«(?)=  f  if/n-i(t')<//i{t-t')dt' ,  (6) 

Jo 

with  The  f„-old  distribution  of  first  sojourn  time 

was  discussed  earlier  in  some  detail  [23].  However,  a  careful 
analysis  of  Eq.  (5)  can  help  the  reader  to  realize  the  rationale 
behind  this  distribution  of  sojourn  times.  The  first  term  on 
the  right-hand  side  of  Eq.  (5)  corresponds  to  the  case  when 


Pi(t)-p2(t) 

Pi(0)~p2(0)' 


(ID 


as  one  can  easily  check  by  Laplace  transforming  both  sides 
and  using  the  GME  with  the  fa-old  memory  kernel  (see  Sec. 
II  B).  old  memory  kernel.  In  other  words,  the  correlation 
function  of  arbitrary  age  mirrors  the  extension  of  the  On- 
sager  principle,  which  is  usually  limited  to  infinitely  aged 
systems  [7],  to  physical  conditions  of  arbitrary  age.  It  has  to 
be  pointed  out  that  at  this  stage  there  is  no  guarantee  that  the 
Onsager  principle  of  arbitrary  age  holds  true.  However,  in 
Sec.  II  B  we  show  that  Eq.  (10)  yields  the  exact  result  of 
Godreche  and  Luck  [14]  for  the  corresponding  correlation 
function.  Furthermore,  in  Sec.  II  C  we  establish  the  same 
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result  using  arguments  based  on  trajectories  rather  than  den¬ 
sities,  thereby  affording  an  independent  construction  of  the 
exact  expression  for  the  fa-old  correlation  function.  All  this 
can  be  thought  of  as  a  compelling  demonstration  of  the  cor¬ 
rectness  of  Eq.  (11),  extending  the  Onsager  principle  to  con¬ 
ditions  of  arbitrary  age. 


B.  Derivation  of  the  exact  expression  proposed  by  Godreche 
and  Luck:  The  probability  perspective 

To  establish  that  the  proposed  approach  yields  the  exact 
expression  of  Godreche  and  Luck  [14],  let  us  express  the 
f„-old  correlation  function  through  the  probability  vector  p 
=  {p\ ,p2)  for  the  dichotomous  variable  £=  +  1  to  have  either 
positive  (state  1)  or  negative  (state  2)  values.  We  have  that 

(t ,  (mm), 

W  = - TAT - =Pl(0)pi(f|U=  0)  +p2{0)p2(t\2,t=  0) 

(s  / 

- P i ( 0)p2(t\  1  ,t  =  0)  - p2(0)pi(t\2,t  =  0) ,  (12) 

where  pj(t\k,t=0)  is  the  conditional  probability  that  the  vari¬ 
able  £  is  in  the  state  j  at  time  f,  given  that  at  time  t=0  it  was 
in  the  state  k.  This  means  that  pj(t\k,t=0)  is  obtained  letting 
those  trajectories  evolve  that  at  time  t=  0  had  £  in  the  state  k. 
For  a  straightforward  evaluation  of  pj(t\k,t=0),  we  use  the 
GME  formalism,  adapted  to  the  f  ,,-old  system,  and  we  take 
into  account  the  initial  condition  pl(Q\k,t=Q)  =  Sik.  Accord¬ 
ing  to  the  GME  the  components  of  the  conditional  probabil¬ 
ity  vector  are  determined  by 

d  ^ 

—Pj(t\k,t=  0)  =  -  dt'®,(t-t')'2,Kjipi(t'\k,t  =  0), 

at  J0  i=l 

(13) 

with  <t>,  (n)  given  by  Eq.  (9)  and  the  elements  of  K  given  by 
Eq.  (3)." 

By  Laplace  transforming  (13)  and  doing  some  algebra,  we 
obtain 


’Mk,0)  =  2  [«I  +  (m)K]7V,.(0M).  (14) 


i=  1 


Defining  the  matrix 


it  +  <f>,  (m)  T>,  (n) 


J  =  [«I  +  <f>,  (h)K]- 


u\u  +  2<l>,  («)]  u\u  +  2<f>,  (m)] 

a  a 

<f>,  (m)  It  +  <f>,  ( u ) 

_ _a _  _ _a _ 

\  u\ii  +  2<J>,  (m)]  u[u  +  2<t>,  (m)]  j 


and  using  the  initial  condition  /?,(0|k,r=0)  =  St  j.,  we  obtain 
for  the  Laplace  transform  of  the  conditional  probability  vec¬ 
tor  the  following  expression: 


Pj(u\k,0)  = 


[u  +  $,  (u)]Sjk  +  <£>,  ( u)(S 

i,k+ 1  +  8j,k- 1) 

a  J  ci  J  J 

u[u  +  2&,  (n)] 


(15) 


Using  Eq.  (15)  we  can  Laplace  transform  (12)  to  obtain 


u  +  <t>,  (it)  u  +  <b,  (u) 

«^g(«) = pM - f - +pi(  o) - f - 

u(u  +  2<E>,  («))  u(u  +  2<E>,  (m)) 

(u)  4>t  (u) 

~  P 1  (0) - - -  p2(  0)  - 


u(u  +  2<f>,  (m))  u(u  +  2d>,  (n)) 

a  a 

(16) 

We  note  that  the  probability  is  normalized,  p1(0)+p2(0)=  1. 
Thus,  it  follows  that 


u  +  <l>,  (m) 


(«) 


u(u  +  2<f>,  (n))  u(u  +  2<f>,  (m)) 

a  a 

x[pi(o)+p2(o)] 

l 

u  +  2<f>,  (m) 


(17) 


confirming  the  correctness  of  the  definition  introduced  in  Eq. 
(10).  Substituting  (9)  into  (17)  we  obtain 


1 


1+2- 


4>t  («) 


1  +  tp(u)  -2  if/,  ( u ) 


1  - 


2if/,(u) 

1  +  if^u)  _ 


(18) 


which  coincides  with  the  results  of  Godreche  and  Luck  [14]. 
Furthermore,  the  ratio  of  the  differences  in  probability  is 
determined  in  Laplace  space  by 

2 

pM-poju)  =  y  (J  1,-  -  T2,)p,(0)  =  1  =  |(,j,  , 

/?i(0)-p2(0)  i=1  Pi(0)-p2(0)  u  +  2  4>,(u)  f 


(19) 


As  pointed  out  earlier,  Eq.  (19)  means  that  one  can  extend 
the  Onsager  principle  from  the  infinitely  aged  systems,  for 
which  Onsager  originally  defined  it,  to  systems  of  any  age. 
In  the  latter  case  the  relaxation  is  proportional  to  the  ffl-old 
correlation  function,  not  to  the  infinitely  old,  or  equilibrium, 
correlation  function.  In  summary,  we  discovered  an  Onsager 
principle  of  arbitrary  age,  at  least  in  the  special  case  of  the 
dichotomous  variables  considered  in  this  paper. 


C.  Derivation  of  the  exact  expression  proposed  by  Godreche 
and  Luck:  The  trajectory  perspective 

It  is  possible  to  again  derive  the  exact  result  of  Eq.  (18) 
from  a  different  perspective,  which  will  allow  us,  in  Sec.  IV, 
to  propose  an  analytic  expression  for  the  f ,,-old  correlation 
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function  as  a  function  of  time.  This  expression,  as  we  shall 
see,  is  not  exact,  but  it  is  shown  numerically  to  be  a  very 
good  approximation  to  the  exact  result. 

The  usual  method  of  connecting  the  correlation  function 
<[' j  to  the  waiting-time  distribution,  within  a  trajectory  per¬ 
spective,  is  to  introduce  a  theoretical  waiting-time  distribu¬ 
tion,  <//(f),  which  cannot  be  observed  directly.  In  fact,  the 
experimental  waiting-time  distribution,  namely,  the  distribu¬ 
tion  of  times  with  alternate  signs,  denoted  by  us  as  £At),  is 
obtained  from  the  theoretical  waiting-time  distribution,  </'  (?), 
by  adopting  the  following  procedure.  We  divide  the  time  axis 
into  bins,  whose  size  is  determined  by  the  waiting-time  dis¬ 
tribution  Then,  these  bins  are  assigned  either  the  value 
1  or  the  value  -1,  by  tossing  a  coin  to  make  the  decision.  It 
is  evident  that  the  intervals  along  the  time  axis  with  the  same 
sign,  are  larger  than  the  time  bins  determined  by  since 
two  or  more  consecutive  coin  tossings  might  have  produced 
the  same  sign.  It  is  shown  [22]  that  the  Laplace  transform  of 
ijt  is  connected  to  the  Laplace  transform  of  i//'  via  the  relation 

nu>  =  — — — ■  (2°) 

2  -  ifj  (u) 

Let  us  use  the  term  event  to  denote  the  coin  tossing  intro¬ 
duced  above.  The  expression  (20)  is  the  result  of  summing 
over  all  possibilities  of  not  changing  sign  with  a  coin  toss, 
which  turns  out  to  be  a  geometrical  series  in  the  Laplace 
representation.  The  correlation  function  (Il  and  the  theoreti¬ 
cal  waiting-time  distribution  function  </f  are  connected 
through  the  relation 


/*co  rt 

dt'ti(t')= i-  dt'fiv). 

S  8  J,  Jo 

(25) 

It  is  straightforward  to  show  that  the  fa-old  experimental 
waiting-time  distribution  and  the  f  ,,-old  theoretical  waiting¬ 
time  distribution  are  connected  through  the  following  sum  of 
convolutions : 

'Ar.M  =  j  S(t)  +  +  “<A  ■  j, 

(26) 

where  the  symbol  ★  denotes  time  convolution.  In  fact,  after  a 
first  interval  of  time  followed  by  a  coin  toss  with  no  change 
of  sign,  determined  by  <//,  /  2,  the  next  intervals  of  time  with 

a 

no  change  of  sign  according  to  the  coin  tossing  prescription, 
are  determined  by  <///2.  The  sum  of  the  convolutions  takes 
into  account  all  the  possible  sequences  of  intervals  of  time 
with  no  change  of  sign  before  a  change  of  sign  of  the  vari¬ 
able  £  eventually  occurs,  and  gives  as  final  result  the  distri¬ 
bution  for  a  first  observed  sojourn  time  t  of  the  variable  £  in 
one  of  its  two  states,  that  is  i/f,  (f). 

Thus  by  summing  the  geometric  series  in  the  Laplace 
variables,  from  (26)  we  obtain 

<A*  («) 

&(«)= - •  (27) 

2  —  ijj  {u) 

Using  Eq.  (20)  we  write  Eq.  (27)  as 


< f>*(0  = 


/•oo 

J,('- 


where  the  average  waiting  time  is  given  by 


/ 

Jc 


(t)  =  Tiji 


(21) 


(22) 


V,  (m)  : 


2<A/a(n) 
1  +  i/Au) 


(28) 


By  Laplace  transforming  Eq.  (25)  and  using  Eq.  (28),  we 
obtain 


&f\u)  =  ■ 


1  -  i }*  (h)  j 


1- 


2i}ta(u) 

1  + 


(29) 


Equation  (21)  determines  that  the  correlation  function  ^(f) 
is  equal  to  the  probability  of  finding  a  window  of  length  t 
without  internal  events. 

The  same  result  can  be  immediately  recovered  using 


namely,  we  again  recover  the  exact  result  of  Godreche  and 
Luck  given  by  Eq.  (18).  This  establishes  the  equivalence  of 
the  trajectory  and  GME  prescriptions  for  this  process. 


ch^)  =  'lC« 


dt'ilUt') 


dt"<P\t"). 


(23) 


We  note  that,  see  Ref.  [7], 

=  7T  f  dt'  i//*(t')  (24) 

ivJ, 

is  the  infinitely  aged  waiting-time  distribution.  Actually,  Eq. 
(23)  is  the  infinitely  aged  correlation  function,  a  special  case 
of  the  more  general  prescription 


D.  Generalized  master  equation  of  arbitrary  age 

We  are  now  in  a  position  to  make  a  preliminary  balance  of 
the  results  obtained  so  far.  The  first  is  that  we  have  general¬ 
ized  the  result  of  an  earlier  paper  [7],  in  that  we  have  derived 
the  GME  of  arbitrary  age 

ypW  =  -  f  &,(t-t')Kp(t')dt',  (30) 

dt  JQ  a 

whose  memory  kernel  <E>,  (t)  is  defined  though  its  Laplace 
transform  by  means  of  Eq.  (9).  We  have  also  shown  that  the 
Onsager  principle,  valid  for  infinitely  aged  systems,  can  be 
extended  to  conditions  of  any  age,  and  that  this  extension 
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allows  us  to  derive  an  exact  expression  for  the  4-old  corre¬ 
lation  function.  However,  the  analytic  results  obtained  so  far 
are  in  the  Laplace  domain.  It  is  desirable  to  achieve  them  in 
the  time  domain,  as  well.  We  now  turn  our  attention  to  the 
latter. 


P-  =  “T>  (33) 

z-  1 

and  the  parameter  T,  characterizing  the  waiting-time  distri¬ 
bution  of  Eq.  (32)  is 


III.  THE  LIOUVILLE-LIKE  APPROACH 


(34) 


In  this  section  we  derive  another  expression  for  the  fa-old 
correlation  function  adopting  a  perspective  where  aging  is 
determined  by  the  out  of  equilibrium  bath  for  the  variable  of 
interest,  i.  There  might  exist  conditions,  as  we  shall  see, 
where  equilibrium  is  not  even  allowed.  The  expression  for 
the  ffl-old  correlation  function  afforded  by  this  perspective  is 
exact,  and  is  thus  equivalent  to  the  Godreche  and  Luck  ex¬ 
pression  of  Eq.  (18).  However,  the  exact  expression  is  im¬ 
plicit,  rather  than  explicit,  and  is  therefore  more  convenient 
for  the  numerical  calculations  done  subsequently. 

Here  we  adopt  the  perspective  of  earlier  work  [13,20]  to 
account  for  the  aging  effects  characterizing  the  fluctuations 
of  the  dichotomous  variable  These  fluctuations  occur 
while  the  environment  of  the  variable,  £,  slowly  drifts.  This 
drifting  process  is  extended  over  time,  and  could  lead  to 
circumstances  where  it  is  not  possible  to  attain  equilibrium 
asymptotically.  In  keeping  with  the  jargon  of  statistical  me¬ 
chanics  the  environmental  or  “irrelevant”  variable  is  called  y 
and  in  the  model  moves  in  the  interval  /=  [0 , 2] .  In  the  semi¬ 
interval  [0,1],  we  use  the  equation  of  motion  for  the  prob¬ 
ability  density 


dp(.y,t) 

dt 


=  -\—yzp{y,t). 
dy 


(31) 


This  is  the  motion  determined  by  a  potential,  with  the  mini¬ 
mum  at  y=  1,  in  the  over-damped  case.  In  the  interval  [1,2], 
the  overdamped  potential  is  the  mirror  image  of  the  potential 
acting  on  the  left  interval.  Consequently,  if  the  initial  condi¬ 
tion  is  located  in  the  internal  part  of  the  interval  [0,1],  the 
particle  moves,  from  the  left  to  the  right,  with  a  deterministic 
motion,  towards  y=  1.  If  the  particle  is  initially  located  in  the 
interior  of  the  interval  [1,2],  it  moves  deterministically  from 
the  right  to  the  left.  When  the  particle  reaches  the  potential 
minimum,  it  is  injected  back,  with  equal  probability,  into  any 
of  the  points  of  the  interval  I  excluding  y=0,  y=  1  and  y  =  2. 
The  time  spent  by  y  within  I  corresponds  to  sojourning  in 
one  of  the  two  states  of  the  variable  £,  either  |+)  or  |— ).  The 
instant  of  back  injection  corresponds  to  the  choice  of  the  new 
state  and,  with  equal  probability,  this  is  either  the  same  state 
or  the  other  state.  The  variable  y  represents  the  environment 
of  the  variable  and  its  initial  distribution  is  given  by  the 
state  of  the  bath.  The  corresponding  waiting-time  distribution 
between  two  consecutive  back  injections  is 

tfr*(t)  =  (ji-  p  (32) 


where  the  index  /jl  is  related  to  z  of  Eq.  (31)  by 


in  accordance  with  the  normalization  constraint.  The  authors 
of  Refs.  [7,18]  have  shown  that  the  essential  properties  of  the 
dichotomous  non-Poisson  fluctuation  can  be  accounted  for 
by  limiting  ourselves  to  this  simplified  picture,  involving 
only  the  semi-interval  [0,1].  In  this  simple  picture  the  aging 
process  is  described  by 

dp(y,t)  d 

J^rL  =  -X--yzp(y,t)  +  \p(ht),  (35) 

dt  dy 

which  takes  into  account  the  back  injection  into  the  semi¬ 
interval,  occurring  with  uniform  probability,  when  the  par¬ 
ticle  reaches  the  point  y  =  1 .  For  simplicity,  but  with  no  loss 
of  generality,  we  fix  \=1. 

Using  the  results  of  Ref.  [20]  the  fa-old  correlation  func¬ 
tion  is  evaluated  as  follows.  The  bath  is  prepared  at  time  t 
=  -ta.  This  means  that  at  time  t--ta  the  distribution  of  y 
within  the  interval  I  is  flat.  In  the  case  z  <  2  this  distribution 
tends  to  the  equilibrium  distribution,  p eq<*  l/yz_l.  If  z>2  this 
distribution  diverges,  thereby  implying  that  the  distribution 
approaches  the  Dirac  delta  function  located  at  y  =  0.  This  is 
the  nonstationary  condition,  the  condition  where  equilibrium 
is  not  allowed,  and  only  a  condition  of  eternal  drift  is  admit¬ 
ted.  Suppose  this  distribution  evolve  for  a  time  ta,  without 
our  observing  it.  This  means  that  we  begin  the  observation 
process  when  the  system  has  a  new  distribution,  different 
from  the  initial  flat  distribution,  and  determined  by  its  time 
evolution  from  t=-ta  to  f=0,  described  by  Eq.  (35).  For  t 
3*0  Eq.  (35)  is  replaced  with  Eq.  (31) 


dp{y,t) 

dt 


—yzp{y,t), 

dy 


(36) 


namely  the  back  injection  process  is  stopped,  thereby  imply¬ 
ing  that  the  population  decreases. 

The  theory  of  Ref.  [20]  relates  the  probability  solution  to 
the  Liouville-like  equation  to  the  fa-old  correlation  function 
as  follows: 


=  f  dypiyj).  (37) 

7o 

Note  that  the  initial  condition  p(y,  0)  is  obtained  from  Eq. 
(35)  moving  from  the  flat  distribution  at  t=-ta.  The  time 
evolution,  corresponding  to  the  observation  process,  is  deter¬ 
mined  by  Eq.  (36).  This  simple  picture  is  a  fair  representa¬ 
tion  of  the  description  made  in  terms  of  trajectories  in  Sec. 
II.  The  fact  that  the  norm  of  p(y,t),  when  p{y  ,t)  is  described 
by  Eq.  (36),  is  not  conserved,  reflects  the  occurrence  of 
jumps  from  one  state  to  the  other,  making  the  population  of  a 
given  state  decrease.  Let  us  remark  that  the  correlation  func¬ 
tion  is  determined  by  the  antisymmetric  part  of  the  whole 
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distribution  [20].  This  formal  condition  corresponds  to  ob¬ 
serving  the  time  evolution  of  p(y,t),  with  y  ranging  from  v 
=  0  to  y=  1,  under  the  action  of  Eq.  (31),  with  no  back  injec¬ 
tion  process.  The  process  of  back  injection  is  essential  to 
determine  the  correlation  function  of  arbitrary  age.  However, 
once  the  parameter  ta  is  fixed,  and  with  it  the  age  of  the 
correlation  function,  the  evaluation  of  the  correlation  func¬ 
tion  is  done  by  imagining  the  bath  frozen  in  the  distribution 
corresponding  to  this  fixed  age.  This  is  the  reason  we  use  Eq. 
(31)  to  determine  the  correlation  function  rather  than  Eq. 
(35). 

The  process  of  back  injection  is  essential  for  the  slow 
environmental  drift,  but  the  correlation  time  of  age  ta  is  de¬ 
termined  by  the  age-fixed  condition,  the  bath  being  ta  old, 
and  keeping  this  age  forever. 

Our  goal  is  to  calculate  the  correlation  function  of  arbi¬ 
trary  age.  According  to  Eq.  (37),  this  goal  requires  that  we 
evaluate  p(y,t),  first.  To  carry  out  this  calculation  we  follow 
the  procedure  illustrated  in  detail  in  Ref.  [24].  The  time  evo¬ 
lution  from  t=-ta  to  t<  0  is  described  by  Eq.  (35),  yielding 
for  the  solution 


p(y,t)  = 


l 


[1  +  (z-l)(t  +  0/“1?/fe-I) 

p(l,r)dT 


f 


[1  +  (z-  l)(f-  r)yz  ']z/(z  1)5 


(38) 


with  p(  1,t)  determined,  in  turn,  by  Eq.  (35). 

As  already  stated,  for  f>0,  p{y ,t)  is  governed  by  Eq. 
(36),  so  that  its  solution  is,  given  the  initial  value. 


p(y,t)  = 


p{y,  0) _ 

[l  +  fc-i)/-1*]*^ 


rj=y!\_  1  +  (z-  l)jz 


(39) 


Now,  we  can  calculate  the  final  expression  for  the  age- 
dependent  correlation  function  given  by  Eq.  (37).  After  some 
simple  algebra  we  obtain  the  following  expression: 


O^f)  = 


1 


[l  +  (z-l  )(t  +  ta)]m~l) 

p{l,r)dT 


£ 


[1  +  (z-l)(?-  T-)]1^-1) 


=  ¥‘(t  +  t, 


j+r 

•J  —i 


p(l,  r)dr 


[l  +  (z-l)(f-r)]1/fe-1) 

=  ^(f+fa)  +  [i-^(0] 

p(l,  r)dr 


i: 

/: 


[l  +  (z-l)(f-r)]1/(z-1) 
p(l,  r)dr 


I -ta  [  1  (z  l)r]1/(z_1) ' 

Note  that  in  Eq.  (40)  we  are  using  the  definition 


(40) 


roo 

’fr’f(/)  =  J  ltf{T)dT.  (41) 

This  notation  for  the  survival  probability  is  also  used  in  Sec. 
IV  A.  We  notice  that,  t >0  and  p{\  ,  f)  is  calculated  from  Eq. 
(35)  with  y=l  and  —ta<t< 0.  This  is  an  exact  expression 
that  turns  out  to  be  convenient  to  numerically  check  the  pre¬ 
scriptions  of  Sec.  IV. 

It  is  worth  using  the  perspective  of  this  section  to  shed 
further  light  into  the  problem  under  discussion  in  this  paper. 
Let  us  notice  the  following. 

If  the  age  vanishes  ffl— >  0,  then  the  f  ,,-old  correlation  func¬ 
tion  reduces  to  the  probability  of  no  event  occuring  <t>l,a)(f) 

If  the  age  increases  without  limit  s-°°,  then  the  fa-old 
correlation  function  reduces  to  a  known  result  cI)</'l*(t) 
where 

^  =  [l  +  (z-l)f](2-z)/(2-b  (42) 

is  the  usual  correlation  function  calculated  at  equilibrium. 
So,  we  must  observe  a  crossover  between  two  distinct  re¬ 
gimes.  Moreover,  we  anticipate  that  we  recover  these  results 
from  a  different  perspective  in  Sec.  IV. 

The  perspective  adopted  in  this  section  makes  it  easier  to 
understand,  on  physical  grounds,  why  the  non-Poisson  con¬ 
dition  produces  aging.  The  Poisson  condition  corresponds  to 
z=l  in  the  nonlinear  dynamical  process  of  Eq.  (35),  and  in 
this  case  the  distribution  is  always  flat.  Aging  is  the  process 
of  slow  regression  of  the  bath  variable  to  equilibrium,  a  pro¬ 
cess  that  can,  in  principle,  last  forever.  With  this  representa¬ 
tion  it  is  possible  to  study  the  correlation  function  of  arbi¬ 
trary  age  also  in  the  case  z>2.  Allegrini  et  al.  [20]  give  a 
detailed  discussion  of  the  effect  of  external  perturbations. 

In  the  case  z<2  equilibrium  is  possible.  An  external  per¬ 
turbation  can  be  used  to  create  a  condition  whereby  the  sys¬ 
tem  is  out  of  equilbrium.  The  observation  of  the  population 
difference  between  the  two  states,  £=  1  and  £=— 1,  is  equiva¬ 
lent  to  determining  the  correlation  function.  Note  that  the 
initial  condition  that  we  use  to  study  the  time  evolution  of 
p(y,t)  during  the  observation  process,  Eq.  (36),  is  the  anti¬ 
symmetric  distribution  of  the  interval  [0,2]  folded  into  [0,1] 
(see  Ref.  [20]).  It  is  evident  that  in  practice  the  equilibrium 
correlation  function  is  never  observed.  In  fact,  the  external 
perturbation  should  create  a  nonsymmetric  initial  condition 
in  the  whole  interval  [0,2].  The  antisymmetric  portion  of  this 
distribution,  folded  into  the  interval  [0,1],  should  create  the 
same  out  of  equilibrium  distribution  as  that  produced  by 
folding  in  the  same  way  the  anti-symmetric  portion  of  the 
equilibrium  distribution.  As  pointed  out  in  Ref.  [20],  this  is 
impossible  to  do  in  practice  with  a  a  perturbation  lasting  for 
a  finite  time,  and,  thus,  with  any  realistic  experimental  pro¬ 
cedure. 

IV.  AN  ANALYTICAL  EXPRESSION  OF  THE 
CORRELATION  FUNCTION  OF  ARBITRARY  AGE 

The  exact  expressions,  Eq.  (18)  and  Eq.  (29),  allow  us  to 
determine  the  dependence  of  the  age-dependent  correlation 
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function  on  time  via  inverse  Laplace  transforms.  This  inver¬ 
sion  is  done  numerically  since  we  are  unable  to  carry  out 
direct  analytic  inversion  of  the  Laplace  expression.  There  is 
therefore  the  need  to  find  an  analytical  expression  for 
<f>£a)(f),  even  if  the  price  we  pay  is  an  approximation  and  the 
loss  of  exactness.  This  section  is  devoted  to  the  derivation  of 
an  accurate,  albeit  approximate,  analytical  expression  for  the 
correlation  function  of  arbitrary  age. 

A.  An  exact  implicit  expression  for  the  correlation  function 
of  arbitrary  age 

We  write  namely  the  probability  of  not  finding 

events  inside  the  interval  (0,f),  given  the  occurrence  of  an 
event  at  time  —ta,  in  the  following  implicit  form: 

®£J(f)  =  +  ta)+\“  dt^(ta  -  .  (43) 

do 


tained  from  0(f):  in  fact,  if  in  Eq.  (5)  we  make  the  substitu¬ 
tion  0(f)— >  0*(f),  we  obtain  exactly  Eq.  (45),  i.e.,  the  defini¬ 
tion  of  0  (t). 

la 

B.  Equivalence  with  the  exact  expression  of  Godreche 
and  Luck 

As  we  said,  Eq.  (43)  is  exact.  Consequently,  this  equation 
is  expected  to  be  equivalent  to  the  exact  expression  proposed 
by  Godreche  and  Luck  in  Ref.  [7].  In  this  subsection,  as  a 
double  check,  we  prove  that  Eq.  (43)  is  in  fact  equivalent  to 
the  expression  of  Godreche  and  Luck,  and  To  do  that,  let  us 
differentiate  both  sides  of  this  equation  with  respect  to  t  to 
obtain 

— <J>£a)(f)  =  -  0*(f  +  ta)  +  f  dt,i//\ta  -  tj)  . 

dt  J  q  dt 

(46) 


This  is  the  probability  of  finding  an  event  at  a  time  larger 
than  t  from  the  preceding  one.  The  right-hand  side  of  Eq. 
(43)  is  the  sum  of  two  terms,  the  first  corresponds  to  no  event 
occurring  in  the  whole  interval  (-fa,f)  and  the  second  term 
takes  into  account  the  possibility  of  at  least  one  event  occur¬ 
ring  between  —ta  and  0.  The  instant  -fz  signals  the  first  of 
these  events  (or  the  event,  if  there  is  only  one).  The  coin¬ 
tossing  procedure  has  the  double  effect  of  rejuvenating  the 
process  (the  correlation  function  now  has  age  tj)  and  of  fac¬ 
toring  the  two  probability  functions  inside  the  integral. 

It  is  interesting  to  notice  that  Eq.  (43)  is  exact,  which  can 
be  verified  by  differentiating  both  sides  of  the  equation,  and 
making  use  of  (25)  and  (41).  This  procedure  yields 

0* (t)  =  0*(f  +  ta)+  f  c/f10i,'(ffl-f1)0*(f),  (44) 

°  do 

namely  an  implicit  (but  analytic)  expression  for  0,  .  A 

la 

straightforward  iteration  of  (44)  leads  to 

0*W  =  0*(f  +  4)+  [  dt\ifj*(ta-  tj)i/r*{t  +  L) 

“  do 

+  [  dtiif/*{ta-tj)  f  dt2if/*(tl-t2)i//*(r  +  t2) 

J  0  Jo 

+  [  dt]i(ii(ta-t])  f  dt2if/*(t,  -  t2) 

Jo  Jo 

xf  dt3if/*(t2-t3)i//*(t  +  t3)+  ■■■  ,  (45) 

do 


A  more  tractable  expression  is  obtained  by  taking  the 
Laplace  transform  of  this  new  expression  with  respect  to  ta  to 
obtain 


Jo 


0(s)-  dye  syijj  (y) 


(47) 


from  which  it  follows  that 


f 

Jo 


es,[0»-  dy  e-'V(>0] 


1  -  i//(s) 


(48) 


The  right-hand  side  of  this  equation  appears  formidable,  but 
one  can  easily  see  by  Laplace  transforming  Eq.  (44)  with 
respect  to  ta,  that  it  is  nothing  more  than  the  Laplace  trans¬ 
form  with  respect  to  ta  of  -0*(f),  namely  -0*(f).  Conse¬ 
quently,  by  inverse  Laplace  transforming  Eq.  (48)  we  obtain 


(49) 


which,  by  integration,  yields 


<D^(f)  =  <D^(0)  - 


(50) 


The  initial  value  of  the  fa-old  correlation  function  can  be 
determined  by  explicitly  Laplace  transforming  Eq.  (43),  cal¬ 
culated  at  t= 0,  with  respect  to  ta,  thereby  yielding 


i.e.,  a  sum  over  all  possibilities  of  finding  events  at  both 
times  1  and  —ta  with  no  events  between  0  and  t.  This  is 
exactly  the  definition  of  0,  (t).  In  detail,  the  first  term  in  the 
right-hand  side  of  (45),  corresponds  to  the  case  in  which  no 
event  occurs  in  the  interval  l=(-ta,0],  while  all  further 
terms  refer,  respectively,  to  one  event  at  f ,  in  /,  to  two  events, 
at  f |  and  t2,  in  /,  and  so  on,  with  tl<t2<t3<-  ■  Note  that 
0*  (t)  is  obtained  from  if/*(t)  in  the  same  way  0,  (f)  is  ob- 

a  a 


<J.6»(0)  =  'V(s)  +  0*(s)<J>^(O),  (51) 

which  simplifies  to 

[i-P(s)]6?(p)  =  4'\s). 

On  the  other  hand,  we  know  that  AL*(.s,)  =  [l  -  0*(v)]/v,  so  the 
initial  Laplace  variable  is 
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4>f(0)=l/s. 

It  is  obvious  that  the  inverse  Laplace  transform  of  this  last 
expression  yields  the  initial  condition  for  the  correlation 
function  (0)=  1,  as  it  must  be,  based  on  the  definition  of 
Eq.  (1)  when  tx  =  t2  and  £=  +  l,  a  dichotomous  variable.  Us¬ 
ing  this  initial  condition  in  (50),  it  follows  that 

= i  -  (*!)=*;«,  (52) 

Jo  a 

which,  as  we  know  from  Sec.  II,  coincides  with  the  exact 
expression  for  the  ffl-old  correlation  function  of  Ref.  [14]. 


C.  Approximation  through  iterative  expansion, 
and  truncation 

Although  we  established  that  Eq.  (43)  is  exact,  due  to  its 
implicit  nature,  it  cannot  be  used  to  obtain  an  analytic  func¬ 
tion  without  carrying  out  a  reasonable  approximation.  The 
approximation  we  select  rests  on  replacing  <I>^(f)  in  the 
right-hand  side  of  Eq.  (43)  with  the  function  A^(t),  which 
denotes  a  correlation  function  of  uncertain  age,  which  is, 
however,  younger  than  the  ffl-old  correlation  function  defined 
by  Eq.  (43).  This  correlation  function  reads 


i~HO  ' 


(53) 


The  easiest  way  to  derive  (53)  is  to  adopt  the  language  of 
conditional  probabilities.  The  numerator  of  this  expression  is 
the  probability  of  not  finding  an  event  between  0  and  f,  mi¬ 
nus  the  probability  of  not  finding  an  event  between  —ta  and  t. 
If  we  call  A  the  condition  of  no  event  in  the  interval 
(0,f),  and  B  the  condition  of  at  least  one  event  in  the  interval 
(—4,0],  then  the  numerator  of  the  right-hand  side  of  (53)  can 
be  identified  with  the  joint  probability  P(A,B).  On  the  other 
hand,  the  denominator  of  (53)  is  simply  P(B),  the  probability 
of  the  event  B  occuring,  and  therefore  the  function  A('"Ht) 
can  be  identified,  as  it  should,  with  the  conditional  probabi- 
lily  P(A\B),  namely  with  the  probability  of  finding  no  event 
between  0  and  t  (i.e.,  it  is  a  correlation  function),  given  an 
event  between  -ta  and  0  (i.e.,  it  has  an  age  younger  than  ta). 

Finally,  by  plugging  Eq.  (53)  into  Eq.  (43)  we  obtain 


=  ¥*(f  +  O  +  [1  - 


(54) 


FIG.  1.  The  ta- old  correlation  function,  O^(r),  for  different  val¬ 
ues  of  ta.  Curves  represent  a  numerical  integration  of  (40)  with  z 
=  5/3  and  T=  1.5,  while  dots  correspond  to  the  approximate  formula 
(54). 

4=o°.  The  correlation  function  dg  ,,(f)  has  a  faster  decay, 
with  decreasing  values  of  ta.  Finally,  the  lowest  curve  repre¬ 
sents  the  correlation  function  with  zero  age,  namely  'E,  (t). 
The  symbols  in  Fig.  1  overlaying  the  continuous  curves  rep¬ 
resent  the  calculations  using  Eq.  (54)  with  various  ages.  We 
see  perfect  agreement  for  ta =£ (t)  and  for  ta>(t).  The  agree¬ 
ment  remains  good  for  intermediate  values,  with  discrepan¬ 
cies  comparable  to  the  numerical  round-off  errors. 


V.  PHYSICAL  CONSEQUENCES  OF  THE  RESULTS 
OBTAINED 

The  results  of  the  earlier  sections  establish  that  the  GME 
of  arbitrary  age,  proposed  in  this  paper,  Eq.  (30),  fits  the 
requirement  of  complete  equivalence  with  the  current  litera¬ 
ture  on  non-Poisson  renewal  processes.  It  is  important  to 
stress  that  a  significant  model,  generating  non-Poisson  statis¬ 
tics  from  a  renewal  process,  has  been  proposed  by  Bouchaud 
and  others  [25,26].  This  model  has  raised  wide  interest  and 
has  recently  been  applied  to  a  variety  of  pheonomena 
[27-31]. 

The  Bouchaud  trap  model  refers  to  the  spin-glass  dynam¬ 
ics,  which  is  modeled  with  a  set  of  branches,  with  different 
energies  E,  having  a  distribution 

p{E)  =  (3g  exp(-  f3sE) ,  (55) 


This  expression  is  not  exact,  but  it  is  analytic.  Moreover,  the 
form  of  Eq.  (43)  suggests  an  iterative  approach,  and  we  can 
therefore  refine  our  result,  by  replacing  cLg(r)  with  A*g(f), 
after  an  arbitrary  number  of  iterations.  In  Fig.  1  we  check  the 
accuracy  of  the  first-order  approximation  to  the  exact  expres¬ 
sion  given  by  (54),  by  comparing  it  to  the  exact  prediction  of 
Eq.  (40).  The  curves  in  Fig.  1  represent  a  numerical  treat¬ 
ment  of  (40)  with  different  values  of  ta.  In  this  example  /jl 
=  2.5  and  T=1.5  (as  said,  A.=  l).  The  highest  curve  in  the 
figure  represents  the  stationary  correlation  function,  namely 


where  Tg=  1  / /3„  is  the  glass  transition  temperature.  For  each 
branch  this  model  assumes  the  Arrhenius  prescription 

t  =  t0  exp(fJE) ,  (56) 

where  T=l//3  denotes  the  sample  temperature,  and  t  the 
time  necessary  to  overcome  the  barrier  of  intensity  E,  by 
means  of  thermal  fluctuations. 

Thus,  using  the  relationship  between  the  energy  distribu¬ 
tion  and  the  waiting-time  distribution. 
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ip(t)dt  -  p{E)dE ,  (57) 

after  some  simple  algebra,  these  authors  [25,26]  obtain  an 
inverse  power-law  waiting  time  distribution 

•Pit)  =  (58) 

where  the  power-law  index  is  given  by 

v =^.  (59) 

P 

For  /3  >  (3,,  the  resulting  waiting-time  distribution  asymptoti¬ 
cally  corresponds  to  that  of  Eq.  (32),  with  /i<  2.  In  this  case, 
the  stationary  GME,  namely,  the  infinitely  aged  master  equa¬ 
tion,  does  not  exist,  and  we  have  to  consider  conditions  for 
the  system  always  being  out  of  equilibrium.  The  approach 
discussed  herein  does  not  have  any  limitation,  even  if  it  is 
inspired  to  original  work  done  moving  from  the  stationary 
condition  [see  Eq.  (7)]. 

It  is  important  to  stress  that  there  exists  increasing  interest 
in  understanding  the  blinking  quantum  dot  phenomena.  In 
addition  to  the  paper  by  Jung  et  al.  [10],  we  mention  here 
more  recent  references  [9,11,32-37].  The  work  of  Brokman 
et  al.  [9]  proves  that  these  processes  undergo  the  same  form 
of  aging  as  the  Bouchaud  trap  model,  thereby  involving,  also 
in  this  case,  the  non-Poisson  renewal  perspective.  It  is  im¬ 
portant  to  stress  that  in  this  field  the  recent  work  of  Veberk  et 
al.  [37],  addressing  explicitly  the  issue  of  the  intermittent 
fluorescence  of  blinking  quantum  dots,  rests  on  a  model  that 
turns  out  to  be  equivalent  to  the  Bouchaud  trap  model.  In 
fact,  according  to  [37],  the  electron  responsible  for  the  fluo¬ 
rescent  emission  is  trapped  at  a  distance  r  with  probability 

p(r)  =  a  exp(-  apr) .  (60) 

The  recovery  probability  is 

k(r)  =  ( l/f0)exp(-  akr) ,  (61) 

implying  that  the  electron  is  trapped  for  a  time  1  given  by 

t=t0exp(akr).  (62) 

It  is  evident  that  this  model  is  equivalent  to  the  Bouchaud 
trap  model.  To  realize  this  fact  it  is  enough  to  replace  the 
energies  E  of  the  Bouchaud  model  with  the  distances  r  of  the 
model  of  Ref.  [37].  The  conclusion  is  that  also  in  this  case  an 
inverse  power  law  is  generated, 

iA0) =  v [TiT’  (63) 


where 


i/  —  .  yv_»-r ) 

ak 

If  <  ak,  we  obtain  a  waiting  time  distribution  with  a  di¬ 
verging  first  moment,  corresponding  to  the  blinking  quantum 
dot  condition  [37]. 

To  move  from  the  physics  of  spin  glasses  to  the  physics  of 
quantum  dots,  as  we  plan  to  do  in  future  publications  with 
the  help  of  the  results  of  this  paper,  it  becomes  important  to 


address  fundamental  issues  such  as  the  emission  and  absorp¬ 
tion  process.  This  leads  us  naturally  to  using  the  GME  we 
derived  in  an  previous  section.  If  an  out  of  equilibrium  con¬ 
dition  is  created  at  a  given  time  ta,  then  after  the  preparation 
of  the  material,  we  predict  the  relaxation  process  of  the  quan¬ 
tity 

U(t)=pl{t)-p2{t)  (65) 

to  be  driven  by 

d  f 

-n(f)  =  -  dt'<b,  (t-t'W)dt'.  (66) 

dt  J0 

It  is  important  to  stress  that  this  result  rests  on  a  GME  that 
has  an  origin  different  from  the  prime  principle  procedure 
adopted,  for  instance,  by  Zwanzig  [12],  being  essentially 
based  on  experimental  observation.  Is  this  compatible  or  not 
with  a  Liouville  or  Liouville-like  treatment?  This  is  an  im¬ 
portant  issue  that  is  discussed  in  the  next  section. 


VI.  CONCLUDING  REMARKS 

The  central  result  of  this  paper  is  the  generalization  of  the 
GME  discussed  in  Ref.  [7]  to  one  of  arbitrary  age.  Another 
result,  closely  related  to  the  GME  of  arbitrary  age,  is  the 
generalization  of  the  Onsager  principle  to  physical  condi¬ 
tions  of  any  age.  The  validity  of  this  generalization  of  the 
Onsager  principle  is  confirmed  by  the  fact  that  Eq.  (19),  the 
generalized  Onsager  principle,  yields  Eq.  (18),  and  this  is 
shown  to  be  equivalent  to  the  exact  prescription  of  Godreche 
and  Luck,  which  is  independently  rederived  in  Sec.  II  C. 
Another  interesting  result  is  given  by  Eq.  (54).  This  is  an 
analytical  expression  for  fa-old  correlation  function,  whose 
accuracy  has  been  established  using  Eq.  (40)  which  is  an¬ 
other  exact  expression  for  the  ffl-old  correlation  function,  de¬ 
termined  by  the  Liouville-like  approach  of  Sec.  III.  In  the 
special  case  of  non-Poisson  processes  where  the  function 
4>\t)  has  the  form  given  by  Eq.  (32),  the  approximate  ex¬ 
pression  Eq.  (54)  turns  out  to  be  very  accurate.  Of  course 
there  might  be  non-Poisson  processes  where  Eq.  (54)  is  not 
as  accurate  as  in  the  case  presented  here.  However,  the  itera¬ 
tive  procedure  discussed  in  Sec.  IV  C,  allows  us  to  determine 
higher-order  corrections,  should  they  be  necessary  for  a  more 
satisfactory  treatment. 

It  is  important  to  understand  why  Eq.  (54)  is  not  exact,  in 
spite  of  the  fact  that  the  approximation  made  for  its  deriva¬ 
tion  seems  to  fit  the  renewal  nature  of  the  process  under 
study,  where  any  jump  resets  memory  to  zero.  This  is  a  con¬ 
sequence  of  the  infinite  memory  generated  by  non-Poisson 
dynamics,  in  spite  of  the  fact  the  random  events  reset  the 
memory  to  zero.  The  evaluation  of  the  correlation  function 
involves  probabilistic  arguments,  and  with  them  the  infinite 
memory  associated  with  the  probabilistic  treatment  of  non- 
Poisson  processes. 

We  illuminate  the  meaning  of  the  correlation  function  of 
arbitrary  age  by  means  of  the  Liouville-like  density  ap¬ 
proach.  The  observation  of  the  process  of  regression  to  equi¬ 
librium  of  the  population  difference  corresponds  to  evaluat¬ 
ing  the  antisymmetric  distribution,  while  leaving  the 
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symmetric  part  of  the  distribution  free  to  evolve.  If  the  dis¬ 
tribution  is  at  equilibrium,  the  symmetric  part  corresponds  to 
the  equilibrium  distribution  and  the  integral  of  the  left  por¬ 
tion  of  the  antisymmetric  part,  without  back  injection,  re¬ 
gresses  to  zero  as  the  corresponding  equilibrium  correlation 
function.  For  any  other  condition,  the  integral  of  the  left 
portion  of  the  antisymmetric  part  regresses  to  zero  with  an 
analytical  expression  depending  on  the  time  at  which  obser¬ 
vation  begins.  The  regression  continues  as  a  function  of  that 
specific  initial  condition  while  the  symmetric  part  keeps 
moving  towards  equilibrium  independently  of  the  population 
difference.  This  explains  why  the  regression  to  equilibrium 
depends  on  the  initial  condition,  of  any  age,  with  no  further 
dependence  on  the  bath  dynamics  that  keeps  drifting  towards 
equilibrium.  This  also  explains  why  an  emission  or  absorp¬ 
tion  spectrum  [11]  is  not  stationary  and  changes  with  time. 
The  resonant  radiation  establishes  a  connection  between  the 
antisymmetric  and  the  symmetric  parts  of  the  distribution, 
thereby  updating  observation  to  the  changing  bath  condi¬ 
tions. 

It  is  worth  ending  this  paper  with  some  further  remarks 
about  these  theoretical  problems.  We  have  built  up  a  GME  of 
arbitrary  age,  using  an  empirical  approach.  Is  it  possible  to 
derive  the  same  GME  by  using  a  Liouville-like  approach?  In 
principle,  we  should  use  the  Liouville-like  picture  of  Sec.  Ill, 
to  derive,  via  contraction  on  the  bath  variables,  the  same 
GME,  of  arbitrary  age,  as  that  of  Sec.  II  D.  However,  it  is 
evident  that  this  effort,  even  if  we  were  successful,  would  be 
of  limited  help,  for  practical  purposes.  Suppose,  for  instance, 
that  we  have  to  study  the  response  of  the  system  to  an  exter¬ 
nal,  time  dependent,  perturbation.  Would  the  GME  of  arbi¬ 
trary,  but  fixed,  age,  be  useful  for  this  purpose?  It  is  evident 
that  it  would  not.  In  fact,  the  external  perturbation  at  times 
different  from  the  age  of  the  system  would  produce  effects 
departing  from  the  more  realistic  approach  resting  on  per¬ 
turbing  trajectories.  In  the  specific  case  of  the  absorption 
spectrum  of  blinking  quantum  dots  [11]  the  authors,  in  fact, 
adopted  this  trajectory  perspective  to  make  a  theoretical  pre¬ 
diction  that  is  incompatible  with  the  perturbation  of  a  GME 
of  fixed  age. 
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There  already  exists  in  the  literature  at  least  the  discussion 
of  one  case  [19]  that  seems  to  be  a  natural  consequence  of 
this  property.  Sokolov,  Blumen,  and  Klafter  [19]  derived  an 
exact  density  equation  to  describe  a  subdiffusion  process. 
This  equation  corresponds  to  brand  new  initial  conditions. 
This  means  a  condition  where  ta= 0.  Then,  these  authors  per¬ 
turbed  this  equation  with  a  time-dependent  field,  and  found 
that  the  theoretical  result  conflicts  with  the  behavior  of  the 
CTRW  under  the  influence  of  the  same  perturbation.  This 
apparent  contradiction  arises  because  the  time  dependent 
perturbation  corresponds  to  additional  observations,  taking 
place  at  different  time  values,  none  of  them  coinciding  with 
the  observation  time,  but  the  perturbation  at  t= 0.  This  is  true, 
whatever  the  observation  time  is,  either  ta=0,  as  in  Ref.  [19], 
or  °°>fa>0,  a  condition  requiring  the  GME  of  this  paper. 
Regardless  of  the  observation  time  that  we  assign  to  the 
GME,  it  is  impossible  to  make  the  GME  prediction  identical 
to  the  CTRW  prediction,  if  we  require  the  perturbation  to 
remain  external  to  the  system.  The  concept  of  perturbation 
itself  turns  out  to  be  inadequate  to  study  non-Poisson  pro¬ 
cesses,  regardless  of  its  intensity.  Thus  the  results  of  the 
present  paper,  in  addition  to  shedding  light  onto  the  implica¬ 
tions  of  Ref.  [19],  imply  a  violation  of  linear  response  theory. 
The  only  possible  way  to  make  the  density  compatible  with 
the  trajectory  picture  is  to  make  the  external  perturbation 
become  a  part  of  the  system  under  study.  This  means  that  we 
have  to  build  up  a  totally  new,  field-dependent,  GME,  along 
the  lines  of  Ref.  [11].  This  sets  a  limit  on  the  applicability  of 
the  GME  of  arbitrary  age  found  in  this  paper.  However,  this 
result  seems  to  support  the  conclusion  that  the  trajectory- 
density  conflict,  revealed  by  Bologna,  Grigolini  and  West 
[13]  might  be  a  consequence  of  the  aging  properties  emerg¬ 
ing  from  non-Poisson  renewal  process. 
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